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Abstract 

The production of prompt photons is measured in the photoproduction regime of electron- 
proton scattering at HERA. The analysis is based on a data sample corresponding to a total 
integrated luminosity of 340 pb _1 collected by the HI experiment. Cross sections are mea- 
sured for photons with transverse momentum and pseudorapidity in the range 6 < Ej, < 
15GeV and —1.0 < rf 1 < 2.4, respectively. Cross sections for events with an additional 
jet are measured as a function of the transverse energy and pseudorapidity of the jet, and 
as a function of the fractional momenta x y and x p carried by the partons entering the hard 
scattering process. The correlation between the photon and the jet is also studied. The re- 
sults are compared with QCD predictions based on the collinear and on the kx factorisation 
approaches. 
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1 Introduction 



Isolated photons emerging from the hard subprocess ep — > cyX, so called prompt photons, are 
a powerful probe of the underlying dynamics, complementary to jets. Production of isolated 
photons with high transverse momentum can be calculated in perturbation theory. High energy 
electron-proton scattering is dominated by so-called photoproduction processes, in which a 
beam lepton emits a quasi-real photon which either interacts directly with the proton (direct 
process) or fluctuates into partons which then participate in the hard scattering process (resolved 
process). In prompt photon production, the direct process is sensitive to the quark content of 
the proton through the Compton scattering of the exchanged photon with a quark (75 — > 75) as 
depicted in figure Hk). The resolved process (qg — > jq) is sensitive to the partonic structure of 
both the photon and the proton. A typical diagram is shown in figure UJd). Figure [TJ:) and[TJl) 
show typical higher order diagrams. 




Figure 1: Examples of diagrams for the production of prompt photons in photoproduction: a) 
direct and b) resolved interaction with a parton from the proton, c) box diagram and d) radiation 
of a photon from an outgoing quark. 

The HI collaboration has previously measured prompt photon cross sections in photopro- 
duction [0Q| and in deep inelastic scattering (DIS) [0. The ZEUS collaboration has also reported 
measurements of prompt photon production j3]-0. Both experiments found that in photopro- 
duction the inclusive prompt photon cross section is underestimated by next-to-leading order 
(NLO) QCD calculations B6HH1, while there is reasonable agreement for events with a prompt 
photon and a jet (photon plus jet). In DIS, a leading order QCD calculation [9] significantly un- 
derestimates the production of isolated photons and of photons plus jets. NLO predictions ifTOl 
are only available for the latter and also underestimate the cross section. 

This paper presents results of a measurement of prompt photons in photoproduction. The 
data used for the measurement were collected with the HI detector in the period from 2004 to 
2007 and correspond to a total integrated luminosity of 340 pb^ 1 . This amounts to an increase 
in statistics by a factor of three compared to the previous measurement [1J. During this data 
taking period HERA collided positrons or electrons:!] of energy E e = 27.6 GeV with protons of 
energy E p = 920 GeV corresponding to a centre-of-mass energy of ^fs = 319 GeV. 

Isolated photons with transverse energy 6 < Ej, < 15 GeV and pseudorapidityLl —1.0 < 
rf 1 < 2.4 are measured in events with the inelasticity y in the range 0.1 < y < 0.7. This 

1 Unless otherwise stated, the term electron refers to both the electron and the positron. 

2 The pseudorapidity is related to the polar angle 8 as 77 = — lntan(6>/2), where 9 is measured with respect to 
the direction of the outgoing proton beam (forward direction). 
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extends the phase space of previous measurements at HERA towards larger pseudorapidities of 
the photon and to smaller event inelasticities. 

The main background is due to photons produced in hadron decays. For its discrimination 
from prompt photons, various shower shape variables are used. Differential cross sections are 
presented as a function of the transverse energy and pseudorapidity of the photon. For the 
photon plus jet sample, differential cross sections are measured as a function of transverse 
energy and pseudorapidity of the photon and the jet and the momentum fractions x 7 and x p 
carried by the participating parton in the photon and the proton, respectively. Azimuthal angle 
and transverse momentum correlations between the photon and the jet are also studied. The 
cross sections are compared to QCD calculations based on collinear factorisation in NLO (61121 
and to calculations based on the kr factorisation approach |fTT|. 



2 Theoretical Predictions 

The calculation by Fontannaz, Guillet and Heinrich (FGH) j6l|7]| based on the collinear factori- 
sation approach includes the leading order direct and resolved processes — > and their 
NLO corrections. Besides the production of a prompt photon in the hard interaction, photons 
may originate from the fragmentation of a high momentum quark or gluon in the final state. 
The fragmentation process, described by a fragmentation function, is included in the calcula- 
tion as well as the direct box diagram as shown in figure [TJ:). The contribution from quark to 
photon fragmentation to the total cross section of isolated photons is at the level of 4%. The 
contribution from the box diagram amounts to about 10% on average. The calculation uses the 
parton density functions (PDFs) CTEQ6L |[T2]| for the proton and AFG04 lfi~3l for the photon. 
The scales for renormalisation and factorisation //p, are chosen to be /jr = fip = E^. The 
NLO corrections to the LO cross section are significant for the inclusive sample. They increase 
the predicted cross section by a factor 1.15 — 1.42, the corrections being largest at low E% and 
large rf< . For the photon plus jet sample the corrections are much smaller and below 10% on 
average. 

The leading order predictions of Lipatov and Zotov (LZ) ifTTll are based on the kr fac- 
torisation approach. The calculation uses the unintegrated quark and gluon densities of the 
photon and the proton using the Kimber-Martin-Ryskin (KMR) prescription 031 with the GRV 
parameterisations for the collinear quark and gluon densities lfi~5l[T6ll . The kx factorisation 
approach is expected to account for the main part of the collinear higher order QCD correc- 
tions |[Tni . Direct and resolved processes are considered in the calculation, but contributions 
from fragmentation and from the box diagram are neglected. 

To ensure isolation of the photon, the total transverse energy within a cone of radius one in 
the pseudorapidity - azimuthal angle plane surrounding the prompt photon, excluding its own 
energy, is required to be below 10% of Ej, in both calculations. This requirement slightly differs 
from the one used in the data analysis as described in section |4~T] 

The theoretical predictions are compared to the data after a correction for multi parton inter- 
actions, for hadronisation effects and for the different definition of the isolation of the photon. 
The total correction factors f corr are determined with the signal MC described below as the 
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ratios of the cross sections on hadron level with multi parton interactions and the data isola- 
tion criteria, to the cross sections on parton level without multi parton interactions and using 
the cone cut for the isolation of the photon. The correction factors are calculated for each bin 
using the event generators PYTHIA [17] and HERWIG [18J which have a different model for 
hadronisation. The arithmetic means of the two correction factors are used, while half of the 
difference between the two models is taken as the error. The correction factors for the total 
inclusive cross section range from 0.84 to 0.99 with an average of 0.9. They are largest for low 
Ej, and in the forward direction, where the photon isolation is most sensitive to hadronisation 
and to multi parton interactions. The uncertainty of the corrections is typically 8%. 

The leading order MC generator PYTHIA 6.2 IfTTll is used in this analysis for the prediction 
of the signal. The simulation of multi parton interactions fT9ll20l is included. The hard partonic 
interaction is calculated in LO QCD and higher order QCD radiation is modelled using initial 
and final state parton showers in the leading log approximation [21 J. The fragmentation into 
hadrons is simulated in PYTHIA by the Lund string model [|22l . The simulated signal contains 
contributions from direct (figure \T^l) and resolved (figure Q})) production of prompt photons in- 
cluding QED radiation. In addition, processes with two hard partons in the final state (figured]!) 
are simulated. The simulations use the parton densities CTEQ6L lfi~2ll for the proton and SASG- 
1D ll23l for the photon. Different parton density functions for the proton (CTEQ5L [|24l and 
MRST04 fl2l) and the photon (GRV El and AFG04) are used to estimate the influence of 
the parton densities on the predicted cross section, which varies by at most 10%, mainly due to 
changes of the proton PDF. The multi parton interactions reduce the total inclusive cross section 
by 6% on average. The uncertainty of the correction for multi parton interactions is estimated 
by changing the default parameter for the effective minimum transverse momentum for multi 
parton interactions in PYTHIA (PARP(81)) from 1.9 GeV to 1.6 GeV and 2.2 GeV, respectively. 

To estimate the uncertainty of the hadronisation correction, the HERWIG lfT8l generator 
is also used to model the prompt photon signal. HERWIG simulates the fragmentation into 
hadrons through the decay of colourless parton clusters. 

Background to the analysis of prompt photons mainly arises from energetic photons from 
the decay of hadrons like n° and r) in photoproduction events, which constitute more than 90% 
of the total background prediction. Direct and resolved photoproduction of di-jet events used to 
study the background is simulated with PYTHIA. 

All generated events are passed through a GEANT [26J based simulation of the HI de- 
tector which takes into account the different data taking periods, and are subject to the same 
reconstruction and analysis chain as the data. 



3 HI Detector 

A detailed description of the HI detector can be found in [|27ll . In the following, only detector 
components relevant to this analysis are briefly discussed. The origin of the HI coordinate 
system is the nominal ep interaction point, with the direction of the proton beam defining the 
positive z-axis (forward direction). Transverse momenta are measured in the x-y plane. Polar 
(9) and azimuthal ((f)) angles are measured with respect to this reference system. 
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In the central region (15° < 9 < 165°) the interaction point is surrounded by the central track- 
ing system (CTD) , which consists of a silicon vertex detector Il28l and drift chambers all op- 
erated within a solenoidal magnetic field of 1.16 T. The forward tracking detector and the 
backward proportional chamber measure tracks of charged particles at smaller (7° <9< 25°) 
and larger (155° <9 < 175°) polar angles than the central tracker, respectively. In each event 
the ep interaction vertex is reconstructed from the charged tracks. In the polar angular region 
(11° < 9 < 169°) an additional track signature is obtained from a set of five cylindrical multi- 
wire proportional chambers (CIP2k) [|29l . 

The liquid argon (LAr) sampling calorimeter ll30l surrounds the tracking chambers. It has 
a polar angle coverage of 4° < 9 < 154° and full azimuthal acceptance. It consists of an inner 
electromagnetic section with lead absorbers and an outer hadronic section with steel absorbers. 
The calorimeter is divided into eight wheels along the beam axis. The electromagnetic and the 
hadronic sections are highly segmented in the transverse and the longitudinal directions. Elec- 
tromagnetic shower energies are measured with a precision of a(E)/E = 12%/yjE/ GeV©l% 
and hadronic energies with a(E)/E = 50%/v/E/GeV© 2%, as determined in test beam ex- 
periments H311I32II . In the backward region (153° < 9 < 178°), particle energies are measured by 
a lead- scintillating fibre spaghetti calorimeter (SpaCal) ll33l . 

The luminosity is determined from the rate of the Bethe-Heitler process ep — > ep'j, measured 
using a photon detector located close to the beam pipe at z = — 103 m. 

The LAr calorimeter provides the trigger [34] for the events in this analysis. The hardware 
trigger is complemented by a software trigger requiring an electromagnetic cluster in the LAr 
calorimeter with a transverse energy Ej, > 5 GeV. The combined trigger efficiency is about 
85% at E^ of 6 GeV rising to above 95% for E^ > 7 GeV. 

4 Experimental Method 

4.1 Event Selection and Reconstruction 

Events are selected with a photon candidate in the LAr calorimeter of transverse energy 6 < 
Ej, < 15 GeV and pseudorapidity —1.0 < rf 1 < 2.4. Photon candidates are defined as compact 
clusters in the electromagnetic section of the LAr calorimeter with no matching signals in the 
CIP2k. The CIP2k veto rejects candidates, if there is a signal in at least two layers of the CIP2k 
close to the expected hit position. In addition, a track veto is applied for 9 > 45°. It rejects 
candidates, if a track in the CTD extrapolated to the LAr calorimeter front face matches the 
electromagnetic cluster with a distance of closest approach to the cluster's barycentre of less 
than 15 cm. Photon candidates are also rejected if they are close to inactive regions between 
calorimeter modules. 

Neutral current (NC) deep-inelastic scattering (DIS) events are suppressed by rejecting 
events with an electron candidate not previously identified as photon candidate. Electron candi- 
dates are defined as compact electromagnetic clusters in the SpaCal or in the LAr calorimeter. 
In the LAr calorimeter the candidates are required to have an associated track with a distance of 
closest approach of less than 12 cm. The electron suppression restricts the sample to NC events 
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where the scattered electron escapes along the beam pipe in the negative z direction. The low 
electron scattering angle of such events corresponds to virtualities of the exchanged photon in 
the range Q 2 < 4 GeV 2 . In photoproduction the inelasticity y is expressed as y — W 2 / s, where 
W is the 7p centre of mass energy. In this analysis y is evaluated as yh = T l (E—p z )/2E e , where 
the sum runs over all measured final state particles with energy E and longitudinal momentum 
p z . The inelasticity is restricted to 0.1 < yh < 0.7. The cut at low yh removes residual beam 
gas background and the higher cut on y h removes background from DIS events including events 
with prompt photons and events where the scattered electron is misidentified as a photon. This 
background is below 1.5% in the final sample and is considered as a systematic uncertainty. 

In order to remove background events from non-ep sources, at least two tracks are required 
in the central tracker, assuring a good reconstruction of the longitudinal event vertex position 
which is required to be within 40 cm around the nominal interaction point. In addition, topologi- 
cal filters and timing vetoes are applied to remove cosmic muons and beam induced background. 

The shape of the photon cluster candidate is used to further reduce the background. The 
transverse^ radius Rt of the photon candidate is defined as the square root of the second central 
transverse moment Rt = \fjh,, where the fc'th central transverse moment of the calorimeter 
cells distribution is given by = (|f — (r)\ k ). Here, r is the transverse projection of a cell 
position and the averages are calculated taking into account the cell energies as weight factors. 
The requirement Rt < 6 cm reduces background from neutral hadrons that decay into multiple 
photons. In most cases such decay photons are merged into one electromagnetic cluster, which 
tends to have a wider transverse spread than that of a single photon. 

For events where a second electromagnetic cluster is found, the invariant mass M 77 of the 
photon candidate cluster, combined with the closest neighbouring electromagnetic cluster with 
an energy above 80 MeV, is reconstructed. Photon candidates from n° decays where the two 
decay photons are reconstructed in separate clusters are rejected requiring M 77 > 300 MeV. 

Tracks and calorimeter energy deposits not previously identified as photon candidate are 
used to form combined cluster-track objects. The photon candidate and the cluster-track objects 
are combined into massless jets using the inclusive k T algorithm [|35l with the separation pa- 
rameter Ro set to 1. Jets are reconstructed in the pseudorapidity range —2.0 < ?7 jct < 3.0 with a 
transverse momentum of E^ > 4 GeV. Due to the harder kinematical cuts for the photon can- 
didate there is always a jet containing the photon candidate called the photon-jet. All other jets 
are classified as hadronic jets. To ensure isolation of the photon, the fraction z = E^/E^!~ iet 
of the transverse energy of the photon-jet carried by the photon candidate has to be larger than 
0.9. Here, EjJ~ }e is the transverse energy of the photon-jet. This isolation requirement largely 
suppresses background from photons produced in the hadron decay cascade. Only events with 
exactly one isolated photon candidate are accepted. 

For the photon plus jet sample, events are selected with a photon candidate and at least one 
hadronic jet with — 1.3 < i] ict < 2.3. If more than one hadronic jet is selected, the one with the 
highest E^ is used. 

Four additional observables are defined for the photon plus jet sample which are sensitive 
to the underlying partonic process: 

3 In the context of the cluster shape analysis the transverse plane is defined as perpendicular to the direction of 
the photon candidate. 
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The estimators x^° and Xp°, which in the LO approximation correspond to the longitu- 
dinal momentum fractions of the partons in the photon and the proton, respectively, are 
defined as 



(2y h E e ) ^ X * (2E P ) ■ 

These definitions H36U37I reduce infrared sensitivity for x 1 — > 1 compared to the conven- 
tional definition of x 7 = (E^e^ 1 + E^ t e~ r,lct )/(2yE e ). The above definitions make use 
of the energy of the photon only, which has a better resolution than the energy of the jet. 
However, x^° and Xp° may become larger than unity. 

Two observables p± and A<& describe the transverse correlation between the photon and 
the jet, A<P is the azimuthal difference between the photon and the jet, and p± is the 
photon momentum component perpendicular to the jet direction in the transverse plane 

I Pt I 

At leading order the prompt photon and the jet are back-to-back and p± equals zero for 
direct processes. The observable A<P is strongly correlated with p±, but is less sensitive 
to the energies of the photon and the jet. 



The yh, Ej., z and rf et distributions of events with an isolated photon candidate are shown 
in figure |2] together with the MC predictions from PYTHIA for the signal and the background. 
The signal (background) prediction is scaled by a factor 1.45 (1.7) on average. The scaling 
factors vary as a function of r\ as suggested by the cross section measurement (section |5]). In 
all distributions the data are described within errors by the scaled MC predictions. At this stage 
of the analysis there is still a significant contribution of background from the decay products of 
neutral mesons. 



4.2 Photon Signal Extraction 

The photon signal is extracted from the sample with photon candidates by means of a shower 
shape analysis based on the method described in [|2l[38l- It uses the following six shower shape 
variables calculated from the measurements of the individual cells composing the cluster: 



• The transverse radius of the cluster, R T . 

• The transverse symmetry, which is the ratio of the spread of the transverse cell distri- 
butions along the two principal axes. Single photon clusters are expected to be more 
symmetric than multi-photon clusters. 

• The transverse kurtosis, defined as Kt = (a^) 2 — 3, with fx 2 and /i 4 the second and 
the fourth moment of the transverse energy distribution. 
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• The first layer fraction, defined as the fraction of the cluster's energy detected in the first 
calorimeter layer. 

• The hot core fraction, being the fraction of the energy of the electromagnetic cluster 
contained in the hot core of the cluster. It is defined as the energy fraction in four to 
twelve contiguous cells in the first two calorimeter layers, depending on the polar angle. 
The cells include the most energetic cell and are chosen to maximise the energy. 

• The hottest cell fraction, which is the fraction of the energy of the electromagnetic cluster 
contained in the cell with the largest energy deposit. 

The distributions of the shower shape variables are shown in figure |3] for the prompt photon 
candidates with the kinematic cuts as defined above. The shaded band shows the systematic 
uncertainty assigned to the description of the shower shapes as described in section 14.41 The 
data are compared with the sum of the background and the signal MC distributions, which 
describe the data within the systematic error. 

In order to discriminate between signal and background, probability density functions p^ 
for the signal and p\ for the background are defined for each of the six shower shape variables 
i. Simulated events for the signal and the background are used to determine p^ and p\ g . The 
photon and background probability densities are taken as the product of the respective shower 
shape densities with the method described in ||39l . For each event a discriminator D is formed. 
It is defined as the photon probability density divided by the sum of the probability densities 
for photons and background. Figure |4] shows an example of the discriminator distribution for 
the range 0.94 < rf 1 < 1.42 and four different bins in E^. The discriminator has in general 
larger values for prompt photons than for the decay photons. The separation power is decreasing 
with increasing Ej,. The sum of the MC predictions describes the data within the systematic 
uncertainty of the shower shapes. 

Additional event samples are used for the determination of systematic errors related to the 
cluster shapes. The first sample, containing Bethe Heitler events, ep — >■ e'jp, consists of events 
with an electron reconstructed in the LAr calorimeter, a photon in the SpaCal and nothing 
else in the detector. The second, complementary sample, in addition containing deeply-virtual 
Compton scattering P31 events, is selected by requiring an electron in the SpaCal, a photon in 
the LAr calorimeter and no other particle in the detector. These independent event selections, 
denoted BH and DVCS respectively, provide a clean sample of electromagnetic clusters at low 
transverse energies in the LAr calorimeter and are used to study the description of the shower 
shapes of the photons. A third sample is used to monitor the description of the shower shapes 
of clusters initiated by the decay of neutral hadrons. This sample, denoted BG, is background 
enhanced by selecting events with the inverted isolation criteria z = E^/E^J' }Ct < 0.9 and no 
cut on the transverse radius of the photon candidate. 

4.3 Cross Section Determination 

A regularised unfolding procedure [|40] - l44l is used to relate distributions y TCC of reconstructed 
variables (input bins) to distributions x tm e of variables on hadron level (output bins), to deter- 
mine the fractions of signal and background and to correct the data for the detector acceptance. 
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The unfolding matrix A relates the two vectors, Ax true = y TCC . Further details on the method 
can be found in in P4l and are summarised in appendix lAl 

The input is binned in three dimensions in the reconstructed quantities E^, rp and D; the 
latter allows the discrimination of signal and background. The output of the unfolding procedure 
contains the number of signal events N sig in E^-rf 1 bins on hadron level and the amount of 
background events in any of the input bins. Additional underflow and overflow bins are defined 
for each output variable. Therefore the unfolding matrix A also includes migrations into or out 
of the phase space of the measurement. It is computed using signal and background PYTHIA 
simulation. For measurements including jet-related variables, both the input and the output is 
additionally binned in some variable u, where u is E^ 1 , rf et , x 1 , x p , A<E> or p±. 

The stability of the unfolding procedure is checked by varying the number of input bins 
and changing the bin boundaries. The results from the unfolding procedure are compared to a 
bin-by-bin correction method. Agreement is seen within errors for most of the analysis bins. 

Cross sections are presented for Q 2 < 1 GeV 2 . The extracted number of signal events 
in each bin N sig is corrected for a small contribution of DIS events at virtualities 1 < Q 2 < 
4 GeV 2 . For this kinematic region the scattered electron has a non-negligible probability to 
escape detection. If such events contain in addition photons at high transverse momentum, their 
signatures are very similar to the signal process. The corresponding correction factor f DIS is 
determined with the PYTHIA signal MC and is found to be above 0.98 for most of the analysis 
bins. The bin-averaged double differential cross section on hadron level is obtained as 

dV N^g ■ fpis 

where C is the luminosity, AE^, (Arp) is the bin width in Ej, {rp) and N sig corresponds to the 
number of signal events in the bin E^-rp '. Single differential cross sections as a function of 
E^ {rp) are then obtained by summing bins of the double differential cross sections in rp (E^), 
taking into account the respective bin widths. The total inclusive cross section is obtained by 
summing the measured double differential cross section over all analysis bins. The differential 
cross sections in bins of some jet-related variable u is obtained by unfolding triple-differential 
cross sections in E^, rp and u, which then are summed over the bins in Ej, and rp ' . For the 
calculation of cross section uncertainties, correlations between bins are taken into account. 

4.4 Systematic Uncertainties 

The following experimental uncertainties are considered: 

• The measured shower shape variables in the DVCS and BH event samples defined in 
section 1431 are compared to MC simulations. The uncertainty on the shower shape sim- 
ulation for the photon is estimated by varying the discriminating variables within the 
limits deduced from the differences between data and simulation. The uncertainty of the 
description of the background composition and the shower shapes of neutral hadrons is 
obtained accordingly by comparing the shower shapes of the BG event sample with the 
background MC from PYTHIA. The resulting variation of the total inclusive cross sec- 
tion is 11%. The uncertainty varies between 10% and 25% for the single differential cross 
sections increasing towards large rp . 
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• An uncertainty of 1% is attributed to the energy of the photon [|45l and an uncertainty 
of 3 mrad to the measurement of the scattering angle O for events with rf < 1.4. For 
rf > 1.4 the uncertainty is 4% on the energy scale and 4 mrad on the scattering angle. 
The resulting error on the total inclusive cross section is ±1.5%. 

• A 2% uncertainty is attributed to the measurement of the hadronic energy ll44l . The 
corresponding uncertainty of the total cross section ±1%. 

• An uncertainty of ±3% is attributed to the determination of the trigger efficiency. 

• The uncertainty on the CIP2k and track veto efficiency results in an error of ±2.5% on 
the total inclusive cross section. 

• Background from DIS events leads to a systematic uncertainty of 1.5%. 

• An uncertainty in the description of the dead material in the simulation is accounted for 
by varying the probability of photon conversion before the calorimeter by ±10%. For 
polar angles 9 < 20° it is varied by ±30% because of more dead material in the forward 
region. This results in a 1% uncertainty for the cross section measurements in the central 
region and 3% in the most forward rf 1 bin. 

• The ratio of resolved to direct photoproduction events in the MC simulation is changed 
within limits deduced from the measured x 7 distribution B4l . leading to ±1% systematic 
error due to a different acceptance. 

• The luminosity measurement has an error of 3.4%. 

The effects of each systematic error on the cross sections are determined by evaluating an 
alternative unfolding matrix A' using the MC prediction made with the corresponding system- 
atic variation applied. The differences to the default unfolding matrix A'-A are used to evaluate 
the contributions to the error matrices of the results using standard error propagation. The final 
error matrix is split into fully correlated and fully uncorrected parts which are listed in tables [2] 
to |7j The systematic uncertainty obtained on the total inclusive cross section is ±13%. The 
largest contribution to this uncertainty arises from the systematic uncertainties attributed to the 
description of the shower shapes. 



5 Results 

The prompt photon cross sections presented below are given for the phase space defined in 
tabled] 

Bin averaged differential cross sections are presented in figures [5] to [9] and in tables |2] to [7] 
For all measurements the total uncertainty is dominated by the systematic errors. The figures 
also show the ratio of the NLO QCD prediction (FGH) (6H3 to the measured cross section 
R = (J fgh I & mens with the uncertainty of the NLO calculation. The factors f corr (see section[2]) 
for the correction of the theoretical calculations for hadronisation, multi parton interactions and 
the definition of the isolation are given in the cross section tables with their errors. 
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HI Prompt Photon Phase Space 


T 1 " 

Inclusive 
cross section 


6 < E% < 15 GeV 
-1.0 < rp < 2.4 
z = E~±/E~p et > 0.9 
Q 2 < 1 GeV 2 
0.1 < y < 0.7 


Jet definition 


E^ > 4.5 GeV 
-1.3 < rf ct < 2.3 



Table 1: Phase space for the measurement of prompt photon cross sections. Kinematics are 
defined in the HI laboratory frame. 



The measured inclusive prompt photon cross section in the phase space defined in table [T|is 

a{ep ->■ e^X) = 66.9 ± 1.8 (stat) ± 8.7 (syst) pb. 

Both calculations predict lower cross sections of 52.1 + 34 pb (FGH) and 56.7lg^ pb (LZ), 
while the MC expectation from PYTHIA is 46.4 pb. Theoretical uncertainties due to missing 
higher orders are estimated by simultaneously varying and fip by a factor of 0.5 to 2.0. In 
addition, the errors on the theoretical predictions include uncertainties due to the error of f corr 
and due to the PDFs. All these error sources are added in quadrature. 

Differential inclusive prompt photon cross sections dcx/d-E^ and Aa/drp are presented in 
table |2] and in figure [5] The results are compared to a QCD calculation based on the collinear 
factorisation in NLO (FGH) Od, to a calculation based on the k T factorisation approach 
(LZ) IfTTl . Both calculations are below the data, most significantly at low Ej,. The LZ cal- 
culation gives a reasonable description of the shape of rf 1 , whereas the FGH calculation is 
significantly below the data for central and backward photons {rf < 0.9). 

Double differential cross sections d 2 cr /dEj,drf are shown in figure [6] and table [3] for all five 
bins in rf '. The rf bins correspond to the wheel structure of the LAr calorimeter. LZ provides 
a reasonable description of the data with the exception of the lowest Ej, bin in the central rf 1 
(0.2 < rf 1 < 0.9) region. The FGH calculation underestimates the cross section in the central 
(0.2 < rf 1 < 0.9) and backward {rf 1 < —0.6) region. Here, it is significantly below the data. 
The prediction from PYTHIA is also shown. It underestimates the measured cross section by 
roughly 45%, most significantly at low Ej,. 

The prompt photon plus jet cross section is 

a(ep -> e7 jet X) = 50.1 ± 1.7 (stat) ± 6.5 (syst) pb. 

It is similar to the inclusive cross section, since the prompt photon recoils most of the time 
against a prominent hadronic jet. The theoretical calculations predict cross sections of 40.6 1. 19 pb 
(FGH) and 45.71 2a pb (LZ). Both are compatible with the measurement within the errors. The 
PYTHIA expectation of 33.9 pb is again too low. 
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Cross sections for the production of a prompt photon plus jet are presented in figure [7] 
and tables |H [5] as a function of the variables Ej,, rp ', E 3 ^ and if et . Both calculations give a 
reasonable description of the E^ and Ej, cross sections but show deficits in the description of 
the rf et shape. Here, the LZ prediction is too high for jets with rf ct < 0.5, and both calculations 
underestimate the rate of events with forward jets. As in the inclusive case, the FGH prediction 
is too low for rp < 0.2. 

Photon plus jet cross section as a function of the estimators x^° and x^° are shown in 
figure [8] and table |4l Both distributions are described by the calculations within errors. 

Cross sections for the two observables describing the transverse correlation between the 
photon and the jet, p±_ and A<P, are shown in figure [9] and tables |6j |7] Both variables are 
expected to be sensitive to higher order gluon emission. The phase-space is divided into two 
parts: one with x^° > 0.8 where the direct interaction of a photon with the proton dominates 
and one with x^° < 0.8, including significant contributions from events with a resolved photon. 
For x^° > 0.8 both predictions underestimate the tails of the distributions suggesting that there 
is more decorrelation in the data than predicted. For x^° < 0.8 the p± distribution is harder 
than for x^° > 0.8, which reflects the increased contributions from events with a resolved 
photon and from photons radiated from quarks in di-jet events. The FGH calculation poorly 
describes the pj_ distribution but gives a reasonable description of the measurement in A<& for 
x^° < 0.8, except for the highest bin in A<P. The regions A<P — » 180° and p± — > are sensitive 
to multiple soft gluon radiation which limits the validity of fixed order calculations [46] . The 
LZ calculation includes multiple gluon radiation in the initial state before the hard subprocess 
and describes A<P > 170° and p± < 2 GeV, but predicts a significantly lower contribution of 
events in the tails of both distributions as compared to the data. 

The present measurement is compared to the published results of HI [1] and ZEUS in 
the restricted phase space 0.2 < y h < 0.7. For the comparison with the inclusive measurement 
of HI the rf 1 range is restricted to —1.0 < rf 1 < 0.9. For the comparison with the ZEUS 
results for isolated photons with a jet, the kinematic range is changed to 7 < Ej, < 15 GeV, 
6 < E 3 ^ < 17 GeV and —1.6 < if et < 2.4. The results of this analysis are found in agreement 
with the previous measurements [|44l . 



6 Conclusions 

The photoproduction of prompt photons is measured in ep collisions at a centre-of-mass energy 
of 319 GeV with the HI detector at HERA using a data sample corresponding to an integrated 
luminosity of 340 pb _1 . Photons with a transverse energy in the range 6 < E^ < 15 GeV and 
with pseudorapidity —1.0 < rf 1 < 2.4 are measured in the kinematic region Q 2 < 1 GeV 2 and 
0.1 < y < 0.7. Compared to previous measurements, the range of rf is significantly extended, 
and the luminosity of the measurement is increased by a factor three. 

Single differential and double differential cross sections are measured. The data are com- 
pared to a QCD calculation based on the collinear factorisation in NLO (FGH) Od, to a QCD 
calculation based on the k T factorisation approach (LZ) IfTTTl . and to the MC prediction from 
PYTHIA. The predicted total cross section is lower than the measurement by around 20%. Both 
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theoretical calculations underestimate the data at low Ej,. While the LZ prediction describes 
the shape of da/d^ 7 reasonably well, the FGH prediction is significantly below the data for 
backward photons (r? 7 < —0.6). PYTHIA underestimates the data by roughly 45%, most sig- 
nificantly at low Ej,. 

Differential cross sections for photon plus jet are measured as a function of the observables 
Ej,, rp, rf et , Ej,, x^° , and x^° . The measured cross sections as a function of the transverse 
energy of the photon and the jet as well as x^° and Xp° are described within errors by the 
calculations. However, neither of the predictions is able to describe the measured shape as a 
function of rf et . 

Correlations in the transverse plane between the jet and the photon are investigated by mea- 
surements of the difference in azimuthal angle A<P and of the photon's momentum perpendic- 
ular to the jet direction, p±. A significant fraction of events shows a topology which is not 
back-to-back. Neither calculation is able to describe the measured correlations in the transverse 
plane. 

Prompt photon cross section in photoproduction are now measured at a precision of about 
10%, with hadronisation corrections known at the level of 5%. The challenge remains to further 
improve the theoretical calculations and arrive at a deeper understanding of the underlying QCD 
dynamics in this interesting channel. 
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A Unfolding procedure 

The photon signal is extracted using an unfolding procedure to relate distributions y mc of re- 
constructed variables to distributions x true of true variables on hadron level, to determine the 
fractions of signal and background and to correct the data for the detector efficiency. The 
unfolding matrix A which reflects the acceptance of the HI detector relates the two vectors, 
A^true = 2/rec- Each matrix element Ay is the probability for an event originating from bin j of 
Xtrue to be measured in bin % of y rcc . The matrix A is computed using the PYTHIA simulation 
for the signal and the background, interfaced to the GEANT simulation of the HI detector. 

A schematic view of the simplified unfolding matrix A is shown in figure \10\ Each row of 
the matrix corresponds to one element of the vector x tIue - The elements of x trU c are: signal, 
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Figure 10: Schematic sketch of the unfolding matrix for 2 x 2 signal bins and 3x2x5 
reconstructed bins. Two generator cuts (z and y) and the background are taken into account. 
The size of the boxes reflects the number of entries in a bin. 



migration and background bins. Each column of the matrix corresponds to one element of 
the vector y mc . The elements of y Tec are: reconstructed bins and side bins. When solving the 
equation for x tm c the number of efficiency corrected signal, migration and background events 
is determined in one step. 

The input y rec is binned in three dimensions in the reconstructed quantities E^ 2,3 , rj l > 2 and 
D. The binning in D is required for the discrimination of signal and background. Figure [TOl 
shows 3x2x5 "Reconstructed Bins". The signal is binned in the hadron-level quantities E T ' 
and r] A ' B . Figure [101 shows 2x2 "Signal" bins in these variables. 

In addition, x true includes "background" bins in E^ 2,3 and 77 1 ' 2 , in parallel to the recon- 
structed quantities. These bins give the amount of background in each reconstructed bin. The 
background is determined in the unfolding together with the signal contribution. 
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The final unfolding matrix A also takes into account migrations into or out of the phase 
space of the measurement. For each cut on hadron level, used to define the measurement phase 
space (tabled)), a migration bin is added, containing events generated outside of the phase space 
but reconstructed in any of the input bins. In figure [TOl two such "Migr." bins are shown. In 
order to minimise possible biases introduced by the signal MC simulation outside the phase 
space, each migration bin is subdivided into Ej, and rp bins (not shown in the figure). 

The amount of migration from outside of the generated phase space is controlled by includ- 
ing "Side" bins on detector level for each of the "Migration" bins on hadron level. A side bin is 
defined as a narrow slice outside the nominal cut value of the reconstructed variable. The side 
bins are also subdivided into Ej, and rf 1 bins. 

Using matrix A the unfolded distribution a? tmo is obtained from the observed distribution 
y Tec by minimising a x 2 function given by 

2 2 2 2 

X =Xa + t Xli 

where 

X\ = 1/2 • (y TCC - Aft rU e) T V" 1 (y re c - Aftrue) 

measures the deviation of Ax truc from the data bins y rec . Here, V = Cov (yi, yj) is the covari- 
ance matrix of the data, initially approximated by the observed statistical errors. In order to 
avoid a known bias of this procedure [47 J, the unfolding is iterated using an updated covariance 
matrix B41 . constructed from the expected statistical uncertainties. For a given regularisation 
parameter r the regularisation term is defined as \\ = (^truo) 2 - The minimum \ 2 can be calcu- 
lated analytically and is found as 

f true = ((A^A) + lrT^V-Vrec. 

The size of the regularisation parameter r is chosen using the L-curve method [48TI5011. 
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HI Inclusive Prompt Photon Cross Sections 


[GeV] 


da/dE^ uncorr. corr. 
[pb/GeV] 


fcorr 


[6.0,7.0] 
[7.0,8.5] 
[8.5,10.0] 
[10.0,15.0] 


27.24 ±1.86 ±3.04 
12.94 ±0.71 ±1.94 
6.74 ±0.65 ±0.95 
2.02 ±0.17 ±0.24 


n OO + 0.02 

0.89+2-2? 

— U.U4 
n OQ + 0.02 
u - yo -0.04 

96 + - 02 


rp 


da/drf* uncorr. corr. 
[pb] 


f corr 


[-1.00,-0.57] 
[-0.57,0.20] 
[0.20,0.94] 
[0.94,1.42] 
[1.42,2.40] 


18.4 ±1.3 ±2.5 
23.9 ±1.5 ±1.5 
27.7 ±1.2 ±2.3 
19.3 ±1.3 ±3.0 
11.0 ±1.0 ±3.5 


n qq + 0.02 

u ' yy -0.04 

0-94 

0-90 ±g5i 

0.87ig£ 
0-84 + °;£ 



Table 2: Bin averaged differential cross sections for inclusive prompt photon production as a 
function of rf and Ej, in the kinematic range specified in tabled] The bin ranges, the differential 
cross section values, the uncorrected and correlated uncertainties, and f corr are listed, where 
f corr denotes the correction factor applied to the theoretical cross sections. It corrects for multi 
parton interactions, hadronisation and the different algorithm used for the isolation cut. 
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HI Inclusive Prompt Photon Cross Sections 



rp 


El 
[GeV] 


d 2 cr/d-Ejd?7 7 uncorr. corr. 
[pb/GeV] 


f corr 


[-1.00,-0.57] 


[6.00,7.00] 
[7.00,8.50] 
[8.50,10.00] 
[10.00,15.00] 


9.24 ±0.93 ±1.06 
3.75 ±0.57 ±0.44 
1.43 ±0.41 ±0.24 
0.27 ±0.10 ±0.07 


n qq + o.02 

u - yy -0.04 

n qq + 0.02 

u - yy -0.02 

i m+o- 01 

i - ui -0.03 


[-0.57,0.20] 


[6.00,7.00] 
[7.00,8.50] 
[8.50,10.00] 
[10.00,15.00] 


9.19 ±1.04 ±0.59 
5.02 ±0.59 ±0.48 
2.29 ±0.46 ±0.19 
0.76 ±0.12 ±0.06 


n qq + O-02 

u - yo -0.04 

n qA + 0.02 

qfi + o o 2 

u - yu -0.03 

1 00 + 002 

i - uu -0.03 


[0.20,0.94] 


[6.00,7.00] 
[7.00,8.50] 
[8.50,10.00] 
[10.00,15.00] 


10.90 ±0.86 ±0.78 
5.15 ±0.43 ±0.52 
3.28 ±0.34 ±0.24 
0.83 ±0.10 ±0.07 


0.8718S 

0.891^3 
n q4 + °- 01 

u - y ^-0.04 


[0.94,1.42] 


[6.00,7.00] 
[7.00,8.50] 
[8.50,10.00] 
[10.00,15.00] 


7.68 ±1.02 ±1.18 
3.31 ±0.51 ±0.51 
2.27 ±0.44 ±0.38 
0.66 ±0.12 ±0.09 


0.8318SJ 
8 c i + - 02 

u -° d -0.05 

o qn + - 03 

u - yu -0.04 

n qc; + 0-03 

u - yd -0.03 


[1.42,2.43] 


[6.00,7.00] 
[7.00,8.50] 
[8.50,10.00] 
[10.00,15.00] 


4.54 ±0.90 ±1.29 
2.12 ±0.41 ±0.72 
0.86 ±0.31 ±0.35 
0.40 ±0.08 ±0.09 


0.79i^ 

0-821^ 

0-881^ 
0-93 ±™ 



Table 3: Bin averaged double differential cross section for inclusive prompt photon production 
in bins of transverse energy and pseudorapidity of the photon. The errors are correlated between 
different EJ% and different rp bins. More details are given in the caption of table El 



21 



HI Prompt Photon plus Jet Cross Sections 


E% 
[GeV] 


da/cLEj. uncorr. corr. 
[pb/GeV] 


f corr 


[6.0,7.0] 
\7 8 51 
[8.5,10.0] 
[10.0,15.0] 

L ' J 


18.53 ±1.61 ±2.21 
9 93 ±1 06 ±1 39 
5.50 ±0.45 ±0.73 
1.68 ±0.23 ±0.16 


0.82 IS" 
87 + ,-°J 

w,u ' - 0.05 

92 + - 03 

u - yz -0.04 

0-95 + °$ 

— U.Uo 


rp 


da/d?7 7 uncorr. corr. 
[pb] 


f corr 


[-1.00,-0.57] 
[_0 57 o 201 

KJ • \_/ 1 « W ■ KJ 

[0.20,0.94] 
[0.94,1.42] 

L ' J 

[1.42,2.40] 


14.79 ±1.25 ±1.70 
18 57 ±1 47 ±1 75 

_l_ K^J • -.J 1 — 1 — _L • -L 1 — 1 — -L « 1 .J 

21.12 ±1.21 ±1.77 
13.88 ±3.15 ±2.12 
7.31 ±2.73 ±1.54 


n 04 + 0.06 

u - y ^-0.07 

90 + S-2t 

"•^ - 0.05 

n 07 + 0.03 
0-86 + 

— U.U4 

0.841^ 


e& 

[GeV] 


do/dE^ 1 uncorr. corr. 
[pb/GeV] 


f corr 


\A 5 6 21 
[6.2,8.0] 
[8.0,10.0] 
[10.0,15.0] 


« «n _i_n 74 _i_i 40 

6.93 ±1.08 ±0.84 
6.15 ±0.78 ±0.65 
1.88 ±0.33 ±0.17 


85 + ' 04 

u - od -0.05 

n qq + o.05 
90 + 03 

u - JU -0.05 

96 + 03 

u - yu -0.04 




da / d77 jct uncorr. corr. 
[pb] 


f corr 


[-1.3,-0.4] 
[-0.4,0.5] 
[0.5,1.4] 
[1.4,2.3] 


7.1 ±0.7 ±1.0 
14.9 ±0.8 ±2.1 
18.7 ±1.1 ±2.7 
15.3 ±1.2 ±2.2 


0.80 1^ 
0.86 

o.9it8:8* 
o.94i8^ 



Table 4: Bin averaged differential cross section for prompt photon plus jet production as a 
function of E%, rp ', E^ and ?7 jct . More details are given in the caption of table El 
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HI Prompt Photon plus Jet Cross Sections 





do/dx™ 


uncorr. 
[pb] 


corr. 


] ^ corr 


[0.0,0.5] 


23.0 


±3.3 


±3.8 


n 70 + 0.04 
U -' 8 -0.08 


[0.5,0.7] 


44.1 


±9.6 


±8.1 


89 + 06 

u - oy -0.08 


[0.7,0.9] 


70.3 


±12.3 


±10.9 


1 24 + - 04 


[0.9,1.1] 


75.9 


±9.8 


±8.4 


n 00 + 0.04 

U -^-0.04 




da/dx L p ° 


uncorr. 
[pb] 


corr. 


f corr 


[0.001,0.010] 


1257 


±91 


±102 




[0.010,0.025] 


1325 


±65 


±150 


0.88±85S 


[0.025,0.040] 


698 


±90 


±134 


0-90±g^ 


[0.040,0.060] 


341 


±47 


±66 


0-911^ 



Table 5: Bin averaged differential cross section for prompt photon plus jet production as a 
function of and Xp°. More details are given in the caption of table [2j 
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HI Prompt Photon plus Jet Cross Sections 


X 7 


P± 
[GeV] 


da/dp± uncorr. corr. 
[pb/GeV] 


fcorr 


[0.8,1.1] 


[0,2] 
[2,4] 
[4,6] 
[6,8] 


7.75 ±0.50 ±0.32 
3.96 ±0.43 ±0.53 
2.16 ±0.55 ±0.56 
0.60 ±0.53 ±0.36 


0.89t°a 
0.87±HJ 

n oi +0.08 
U-S-L-O.OS 


[0.0,0.8] 


[0,2] 
[2,4] 
[4,6] 
[6,8] 


7.14 ±0.48 ±1.04 
4.52 ±0.49 ±0.80 
2.91 ±0.40 ±0.56 
2.35 ±0.44 ±0.37 


0.84 lg;°? 

u - yi -0.08 

Q6 + 03 

u - yu -0.04 

1 n7 + 0.09 
L -v ' -0.10 



Table 6: Bin averaged differential cross sections for prompt photon plus jet production as a 
function of p±, the photon's momentum perpendicular to the jet direction in the transverse 
plane, separated into two regions with > 0.8 and < 0.8. More details are given in the 
caption of tabled 
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HI Prompt Photon plus Jet Cross Sections 


X 7 


A<P 


da /dA<P uncorr. corr. 
[pb] 


fcorr 


[0.8,1.1] 


[130,150] 
[150,165] 
[165,172] 
[172,180] 


0.19 ±0.04 ±0.01 
0.46 ±0.06 ±0.06 
0.89 ±0.14 ±0.07 
1.38 ±0.04 ±0.17 


0.821°;°* 

0.931°,;°° 
89 + ' 05 


[0.0,0.8] 


[130,150] 
[150,165] 
[165,172] 
[172,180] 


0.27 ±0.03 ±0.04 
0.52 ±0.06 ±0.08 
0.91 ±0.14 ±0.18 
1.21 ±0.11 ±0.15 


u - y ^-0.05 

n qi +0.04 

u - y± -0.08 

0.83l°;° 9 7 
0.84l°;°« 



Table 7: Bin averaged differential cross sections for prompt photon plus jet production as a 
function of A<E>, the difference in azimuthal angle between the photon and the jet, separated into 
two regions with x^° > 0.8 and < 0.8. More details are given in the caption of tabled 
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Figure 2: Distributions of a) the inelasticity y^, b) the transverse energy Ej, of the photon, c) 
the isolation parameter z and d) rf et for events with prompt photon candidates. Data are shown 
as points with error bars indicating the statistical error. The signal prediction from PYTHIA for 
prompt photons is shown as dark full line, and the contribution of direct interactions as a dotted 
line. The background as simulated with PYTHIA is shown as a dashed line. The signal and 
background contributions are scaled on average by a factor 1.5 and 1.7, respectively. The sum 
of the scaled signal and background is shown as the light grey histogram. The vertical dashed 
lines indicate the kinematic region of the cross section measurement. 
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Transverse Radius [cm] Transverse Symmetry Transverse Kurtosis 




First Layer Fraction Hot Core Fraction Hottest Cell Fraction 



Figure 3: Distributions of the shower shape variables that are used to define the discriminant 
for isolated photon identification: a) transverse radius, b) transverse symmetry, c) transverse 
kurtosis, d) the fraction of energy in the first layer, e) the fraction of energy in the hot core 
of the cluster and f) the fraction of energy in the hottest cell of the cluster for all the photon 
candidates. Data are shown as points with error bars for the statistical error. The expectation 
from PYTHIA for the signal (background) is shown as a full (dashed) line. The signal and 
background contributions are scaled on average by a factor 1.5 and 1.7, respectively. The shaded 
band shows the sum of the MC predictions. The width of the band corresponds to the systematic 
uncertainty assigned to the description of the shower shapes. 
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H1 10<E^<15GeV 

• H1 Data 

Sum 

— Signal 
Background 



+ 



0.2 0.4 0.6 0.8 



l"T ■ ". r'-vvi-v'-vi...-,-,.!", 1 -,-', 
0.2 0.4 0.6 0.8 



Figure 4: The distribution of the discriminant D for 0.94 < rp < 1.42 in four different Ej, 
bins. D is used in the identification of prompt photon candidates for events that have passed 
the event selection. Data are shown as points with error bars for the statistical error. The 
expectation from PYTHIA for the signal (background) is shown as a full (dashed) line. The 
signal and background contributions are scaled by a factor 1.5 and 1.7 on average, respectively. 
The shaded band shows the sum of the MC predictions with the systematic uncertainty from the 
description of the shower shapes. 
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Inclusive Prompt Photon Cross Sections 
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Figure 5: Inclusive differential prompt photon cross sections a) da/dE^ and b) dcr/dTy 7 in the 
kinematic range specified in tabled] The inner error bars on the data points indicate the uncor- 
rected error including the statistical error, the full error bars contain in addition the correlated 
errors added in quadrature. The data are compared to a QCD calculation based on the collinear 
factorisation in NLO (FGH) [[6l|71 and to a QCD calculation based on the kr factorisation ap- 
proach (LZ) rTTTl . For each plot the lower figure shows the ratio of NLO QCD to the measured 
cross section (R = ctfgh/c") as a the hatched band. The width of this band shows the un- 
certainty from the NLO calculation only. The data points are shown at R = 1 and their bars 
indicate the experimental uncorrelated uncertainty. The correlated experimental uncertainty of 
the data is indicated by the shaded area. 
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Inclusive Prompt Photon Cross Sections 
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Figure 6: Inclusive double differential prompt photon cross sections d 2 'a / 'dE^drj 1 for five bins 
in rf 1 . The errors between different rf bins are correlated. The kinematic range is specified in 
tabled] The cross sections are compared to theoretical calculations (see caption figure [5]). In 
addition the prediction from PYTHIA IfTTTl is shown as dashed line. 
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Prompt Photon plus Jet Cross Sections 
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Figure7: Differential prompt photon plus jet cross sections a) da/d-Ej, b) da/drj 1 ', c) da/dE; 
and d) da/drf et in the kinematic range specified in tabled] The cross sections are compared to 
theoretical calculations (see caption figure |5]). 
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Prompt Photon plus Jet Cross Section 



-S 100 

Q. 



X 

■a 



80 



■g 60 



40 
20 

1.5 



0.5 r 



H1 



• H1 Data 

y/// < FGH (NLO) 
LZ(k fact) 



'//. 



6444444i j 64444444 



_i I i_ 



a) 




y //////////////'//////, 
corr. error 



* 

;,////////. 



0.5 



.Q 
Q. 

"1500 

-I CL 

X 

■a 

■§1000 



500 



1.5 
1 

0.5 




>■/////'">""'///, 



////'////////////. 



0.02 



0.04 



rLO 



0.06 



Figure 8: Differential prompt photon plus jet cross sections a) da/dx^° and b) da/dXp°, in the 
kinematic range specified in tableQ] The cross sections are compared to theoretical calculations 
(see caption figure |5]). 
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Prompt Photon plus Jet Cross Sections 
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Figure 9: Differential prompt photon plus jet cross sections da/dA<P and da/dpj_ the photon 
momentum transverse to the jet direction, the difference in azimuthal angle between the photon 
and the jet. The kinematic range is specified in table [Q Figures a) and b) show the cross 
section for > 0.8, c) and d) for < 0.8. The cross sections are compared to theoretical 
calculations (see caption figure [3]). 
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